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Abstract 

We present a formalism to calculate the probability distribution function of a 
scalar field coarse-grained over some spatial scales with a Gaussian filter at finite 
temperature. As an application, we investigate the role of subcritical fluctuations 
in the electroweak phase transition in the minimal standard model. It is concluded 
that the universe was in a mixed state of true and false vacua already at the critical 
temperature. 



1 Introduction 



Quantum field theories with spontaneous symmetry breaking usually exhibit sym- 
metry restoration at high temperature and hence our universe is supposed to have 
experienced several phase transitions in the early stage of its evolution as the cosmic tem- 
perature decreased with expansion. Much work has been done on these phase transitions 
and their generic properties such as the transition order and the critical temperature, 
Tc, have been extracted from the effective potential. The dynamics of phase transitions, 
on the other hand, is less understood, partly because it is difficult to clarify the role of 
thermal fluctuations. 

Thermal fluctuations often play a conspicuous role on the dynamics of phase transi- 
tions. For example, at formation of topological defects associated with some symmetry 
breaking 0, magnitude of thermal fluctuations fixes the temperature when the phase 
distributed randomly on each coherent volume becomes stabilized. Also, their statistical 
property is important to set the initial distribution of topological defects. Furthermore, 
thermal fluctuations called subcritical bubbles have been considered important in the con- 
text of the electroweak phase transition, where the one-loop improved finite-temperature 
effective potential shows that it is of first order Gleiser, Kolb, and Watkins 0, 

however, suggested that the barrier between the two minima of the effective potential at 
Tc was so shallow that thermal fluctuations called subcritical bubbles might be dominant 
and that the standard bubble nucleation picture of the first order phase transition might 
not work. 

Thus, it is indispensable to understand the property of thermal fluctuations in order 
to elucidate the real dynamics of the phase transition. Among the property of thermal 
fluctuations, their typical scale and amplitude are particularly important. The former 
gives us the scale we should focus on, and the latter determines whether thermal fluctu- 
ations dominate the dynamics. First, we should notice that in evaluating the magnitude 
of thermal fluctuations, we should pay attention to fluctuations coarse-grained over some 
characteristic coherent scale of thermal fluctuations rather than raw fluctuations. If we 
estimate the typical amplitude out of raw fluctuations, random infinitesimal ones make a 
dominant contribution to them and we cannot obtain a meaningful value to understand 
the role of thermal fluctuations. 

In this paper we present a formalism to calculate the probability distribution function 
(PDF) of a coarse-grained scalar field at finite temperature. The PDF of coarse-grained 
fluctuation has also been studied by Hindmarsh and Rivers ^ and by Bettencourt |0, 
but they gave only the upper limit of probability. 

As an application, we apply our formalism to the electroweak phase transition which 
may be cosmologically important for the possibility of baryogenesis during it [1^. For 



successful baryogenesis, it should be of first order to realize an out-of-equilibrium condition 
||TT| . However, if subcritical fluctuations dominate the dynamics of the phase transition 
it proceeds in a similar manner to a second-order transition and baryogenesis becomes 
impossible. Here we evaluate PDF of the Higgs field on the relevant scales in the minimal 
standard model and investigate whether the universe is homogeneously in the symmetric 
phase at the critical temperature. We have also studied this issue in the previous paper 
pO| , solving a phenomenological Langevin equation for the classical expectation value of 
the Higgs field numerically. But the PDF we formulate here is based on finite-temperature 
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field theory and quantum effects are taken into account more appropriately. 

In the next section, we present a formalism of evaluating the PDF for the coarse- 
grained fluctuations. In Sec. ^ the method established in Sec. || is applied to the elec- 
troweak phase transition. 

2 Formulation of the probability distribution func- 
tion 

We start with the partition function Z of a scalar field (p{x) at the finite temperature 
T = i3~^, defined by 

Z = tip = J Vip{x)p,3[ip{x),(p{x)] , (1) 

where pi3[ip{x),(p{x)] is the diagonal component of the density matrix, pp, which gives 
the relative probability distribution of the configuration ip{x). Then the PDF of the field, 
P[(p{x)], is expressed by path integral as 

= \ V<P{x) exp j d'xCEim]) , (2) 
where x = (x, r), is the Euclidean Lagrangian density given by 

CE[<P{x)] = \{d(t^f + V[<P], (3) 

and the boundary condition Bi implies 0(a3, r = 0) = (t){x, t = f3) = f{x). 

Now we consider the PDF of the field coarse-grained over an arbitrary scale R as, 

^R = J d^xip{x)I{x;R) , (4) 
where I{x] R) is a Gaussian window function with the width R and given by 

Then the probability that ipR is equal to P[(Pr = <^], is formally given by 

P[^n = ^] = J V^{x)P[ip{x)]S{^n - ^) • (6) 
Inserting (H) into this equation, we get the following expression, 

P[ipR = ip] = Vip{x) V<P{x) exp [- j d''xCE[<p{x)]^ 5{ipR - (^) 

^\\b. '^'^^^^ S exp(-za(^) exp [- J ^^^^[^(x)] + p J d^xJ{x)(f){x, 0)) , (7) 
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where J{x) = ial{x; R)/ (3 and the boundary condition B2 imphes 0(a3, r = 0) = 0(a;, t = (3). 
Now, we represent the field, 0(a;,r), in terms of the Matsubara frequency, 

= Lpo{x) + Lphix,r) . (8) 

Here, ipn represents a heavy mode with an effective mass-squared = V^"[v5o] + (^)^ 
and (fih is the collection of heavy modes. In the high temperature regime, the heavy modes 
have large effective masses so that fluctuations of the heavy modes are much smaller than 
that of the zero mode, (po{x), and can safely be integrated out. For this purpose, the 
Lagrangian is expanded around the zero-mode. 



where CE,int{fh', fo) represents higher order terms in (fh which generate multi-loop effects 



and will be omitted hereafter. Here. 



6^C 



'0 

6^Ce 



^ „ is given by 



and for each n mode it acts as 



2rm 



Thus the exponent of integrand of the path integral in (0) becomes 



(10) 



fill 



d'^xCf 



{x)]+p I d^xJ{x)(f){x,0) 
6C, 



(12) 



J d^X I -CeIMx)] -^Y.V*n' 



X 



J{x) 



where 



d^x = I dr I d^x . 



Note that heavy modes also attach to the external fleld J{x). If we flrst evaluate only 
the probability distribution functional P[(f{x)] for the raw fleld ip{x) as was done by 
Hindmarsh and Rivers [§, the above term is omitted in integrating out heavy modes. 
Therefore our formalism is more rigorous than theirs to this end. Also, as a result, the 
propagator in evaluating the connected generating functional deflned below, is the 
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finite-temperature propagator in Hindmarsh and Rivers 0, while in our formalism it is 
the zero-temperature propagator. 

After integrating over heavy modes, the PDF is given up to one-loop order by 



1 /■ r doi 

P[ipR = 0\ = -J V^o{x) J — exp(-ia^) x 



exp 



J{x) + (3 J dxJ{x)(po{x) 



(13) 



where -^[v^o] is the one-loop finite-temperature three-dimensional effective acizon obtained 
by integrating over heavy modes only. In the case ip is coupled to other quantum fields, 
integration over these fields in all modes should also be performed at this stage to obtain 
the effective action of (po- 

Here we comment on the difference between the effective action integrated over only 
heavy modes and that over all modes. We concentrate on the homogeneous part, namely 
the effective potential for simplicity. The effective potential integrated over all modes is 
written by the zeta-function as [ |T^ 



-^allbc 



1 dCois) 



2(3U ds 



s=0 



where U is the box size of the system, or, / d^x, and Co{.s) is defined by 

Co{s) 



n n 



(14) 



(15) 



with a„ a series of positive real discrete eigenvalues of an operator O. In our case, the 

"-E 



operator O is identified with 



Then the corresponding zeta-function becomes 



Co{s) 



r(s) (47r/32)3/2 



rdaa'-'/' y e 



(47r2n2+y"[^o])< 



(16) 



where /i is the renormalization scale and n represents each mode. Since the zero-mode 
contribution to (o{s) is given by 



Co(s)n=0 



2s 



r(g-3/2) 



(17) 



that to — ; becomes 

ds 



dCoiO) 



ds 



n=0 



Gtt 



(18) 



After all, the difference of the two quantities, namely the zero-mode contribution to the 
effective potential, is given by 
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(19) 



We now return to the evaluation of eq. ([131). Below we omit the suffix 0. We ffist note 
the convergency of a-integration in eq. (|T3p. Using the high temperature limit, we find 



E 
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(20) 



Then a-integration becomes 



/ — exp —laif^ — / a x J[x) 
J — oo 27r 2 J ,r, 



/3 



J{x) + p J d^xJ{x)Lpo{x) 



°° da 

— exp 

-00 27r 

°° da 

— exp 



00 



277 



P f 

-ia{^ - ^r) - —a^ J d^xl{x; i?)^ 



-ia{ip - lpr) 



1927r3/2i?3 



(21) 



From the above expression, it is evident that, while a runs from —00 to +00, the significant 
contribution comes from a finite region 



967r3/2i?3 



< a < 



ara{R) ■ 



[22) 



But we do not perform a-integration for the moment. Instead we ffist define the generating 
functional for the connected Green functional W[J] as 



W[J] =/3-Mn|^y" P(^exp + /3 j d^xJ{x)<^{ 

Then the probability distribution function is expressed as 



(23) 



P[^R = ^] = / ^ exp (/3W[J] - + f / d'x Y: J{x) i 



-1 



J{x) 



(24) 



Thus we can obtain the probability distribution P[fR = by evaluating W[J] in princi- 
ple. The generating functional (|23D is often expressed perturbatively since the higher-point 
connected Green function typically involves higher-powers of (small) coupling constants. 
In our case with J{x) = ial{x; R)/P we can formally write 



f3W[J] = [ d^^i---d^XnG^^\x,...Xn)J{xi)...J{x^) 

^ 71/. J 



E 



na 



(25) 



where ((v^ij)") is defined by 
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((<^k)") = J d'xi . . . d'xnG^:^\x^ . . . Xn)I{xi; R)... I{x^; R) . (26) 

As mentioned above, since the dominant contribution comes from a finite range of |a|, in 
many cases it is suitable to use only lower-order terms of (P3|) in (El). 



3 Application to the Electroweak Phase Transition 

Having established a generic formalism, we now apply our formalism to the electroweak 
phase transition. The one-loop finite-temperature effective potential of the Higgs field in 
the minimal standard model is well approximated by 

Vew[<I>] = D{T^ - T^^' - ETcf' + ^A^^^ , (27) 

where D = 0.169 and E = 0.00965. For Mfj = 60GeV, which is the lower bound allowed 
experimentally [^, we find At = 0.035, T2 = 92.65GeV, and the critical temperature 
is given by Tc = 93.39GeV. The temperature dependence of the potential is depicted in 
Fig. At the critical temperature, the effective potential has two minima, at = and 
(p = 0+(Tc) = = 51.5GeV, and the inflection point is equal to 0inf = 10.9GeV. 

Although the shape of the effective potential exhibits a typical feature of a first-order 
phase transition, the potential barrier, which gets smaller for larger M^, is so shallow 
that it has been suggested that the universe was in a mixed state of true and false vacua 
already at T = Tc due to subcritical fluctuations ||^ ^ []T3[-|]TB[ as mentioned in §|I|. 



Furthermore analytical cross-over is observed for Mh > 80GeV in recent results of lattice 
Monte Carlo simulations [|l7| On the other hand, by estimating the amplitude of 

fluctuations on the so-called correlation length, namely the curvature scale at the origin 
of the potential, or on even larger scales, it has also been claimed subcritical fluctuations 



is small enough and the conventional picture of the phase transition works |p |T^ . 

The amplitude of thermal fluctuations depends largely on the coarse-graining scale. 
In our recent work by solving the simple Langevin equation, we have shown that 



although the curvature scale of the potential or the Compton wavelength is a good measure 
of correlation length for a non-selfinteracting massive scalar field in thermal equilibrium, 
it is not sensible to take it as the coarse-graining scale for a field with a more complicated 
potential such as the electroweak Higgs field. We need to determine it from a fundamental 
point of view. 



In ||2^ we have re-considered derivation of the Langevin equation for the expectation 



value of a scalar field at finite temperature originally developed in |21] and [22] using 
non-equilibrium quantum field theory. We have discussed that properties of stochastic 
noise force in the Langevin equation have distinct features depending on whether it arises 
from interactions with bose fields or from fermi fields, and shown that the latter is more 
effective to disturb the system from a homogeneous configuration. Hence we expect the 
correlation properties of these fermionic noises play an important role in determining the 
correlation length of the scalar field. 

As shown in pO|, the spatial correlation function of the stochastic thermal noise force, 
^{x,t), generated by a massive fermion ip through Yukawa coupling, /(pipip, is given by 



6 



n=l 



n=l 



r = \x — y \ ,(28) 



where / is a coupling constant and is the mass of the coupled fermion ip. It damps 
exponentially above the inverse mass scale for m^P ^ 1, and above the scale (3 for 
m^P < 1. In the case with = 0, we find 



TT 



+ 



71-2 cosh i^^TT 



Pr sinh (^tt) sinh^ (^vr 



4^4 ^2 



for r > — . 

TT 



(29) 
(30) 



So the typical damping scale is given by /3/(27r). In both cases the temporal correlation 
damps exponentially if the temporal difference becomes larger than P/{2tt). 

Therefore we should primarily calculate the PDF taking the coarse-graining scale equal 
to this correlation length of the stochastic noise field acting on because it is the only 
source of inhomogeneous evolution of the field. In practice, however, we calculate how 
the PDF changes as we increase the coarse-graining scale from /5/(27r) to the curvature 
scale at the origin of the potential to compare with other work. Note that it is a good 
approximation to treat the noise field as arising from a massless fermion because, even 
if the two phases = and 0+ are completely mixing the average mass of the top 
quark, which is the most strongly-coupled fermion to the Higgs field and provides the 
dominant source of thermal noise, is only about 18 GeV then, much smaller than the 
critical temperature. 

In order to obtain the proper PDF, we need evaluate the generating functional VTfJ]. 
For this purpose, we must calculate the finite-temperature three dimensional effective 
action, which is a formidable task, and so we approximate it by the standard kinetic 
term and the effective potential. That is, the one-loop three-dimensional effective action 
is replaced by 



H\ 



(31) 



Strictly speaking, this replacement is justified only if deviation from homogeneous config- 
uration is sufficiently small. When the field is in a mixed state between the two minima, 
the effective potential should be modified. It is well-known, however, that this modifica- 



tion reduces the potential barrier ||2^, which induces more phase-mixing. Therefore, even 
if we evaluate the probability distribution function by using the above approximation and 
reach the conclusion that two phases are mixed, the conclusion is consistent. 

Below we proceed the discussion by using the above approximation and adopt Vev^fv?] 
as the effective potential Vesf'/']- The alert reader may notice that Kff[</5] should be that 
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integrated only over heavy modes for the Higgs self- interaction. But, since the effect of the 
Higgs self-interaction is so small that it has already been omitted in evaluating Vevi/ 
it makes no change. We fix T = in order to investigate whether the universe is in a 
homogeneous state of the false vacuum at the onset of the phase transition. We also set 
the Higgs mass Mh = 60GeV. Then Vevi^[v?] is given by 

VewM = \m^V^ - ^/^^c^' + ^A^^ , (32) 

where M ^ 0.073Tc and the dimensionless coupling constants read yU ^ 0.058 and 
A ^ 0.21. Since it is impossible to evaluate the generating functional W[J] exactly 
even in the case ( pl|) is adopted, we expand W[J] perturbatively. Also, we neglect the 
gauge-nonsinglet nature of the Higgs field for simplicity and consider only its real-neutral 
component to treat ip as if it was a real singlet field. 

Up to the first order of the coupling constant, W[J] is made up of the graphs as 
depicted in Fig.^. That is. 



+ 21^3 [3 i^xx^xaJa) xa + {^xa^xh^xcJ aJbJ c) xabc\ 

A 



3 \^xx) ^ + 6 {^xx^xa^xhJ a-h) xab + {^xa^xb^xc^xdJ aJhJ cJ d) ^abcd 



4!/33 

-(- (higher — order terms) , (33) 
where is given by 

' J (27r)3fc2 + M2 ^ ' 

Here (• • ■)xab--- iniplies integration over d^xd^ad^b ■ ■ -, and (3 = l/T^. Each loop graph can 
be renormalized into the definition of the normalization of the probability, and the 
mass, M. Therefore, we have only to evaluate three tree graphs. 

First we consider the two-point function. Adding the last term of the exponent in 
eq.(p3), it is given by 



E 



Jd'xJd'yI{x;R)Wj 



d^k 



P^J (27r)3fc^ + M2 + (27™//5): 



I{x;R) 

■ e^<^-yh{x-R) 



27r2 



d^x / d^y 



dk 



d^k 
(2^ 



1 1 

+ 



1 1 

2 ^ e^^ - 1 



-e'^<''-yh{x-R)I{y-R) 



2 e/^"^ - 1 



-I{k;Ry 



(35) 



where /(fc; R) is the Fourier transform of the window function. 
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) 



(36) 



Similarly, the three- and the four-point functions are given, respectively, by 




(37) 



where 



fik) 




(38) 



In order to confirm that it is sufficient to consider graphs with lowest-order in coupling 
constants for the case with the fundamental coarse-graining scale R = /3/(27r), let us 
estimate the magnitude of graphs with n-th order of coupling constants on scale R = 
N(3/{2t[). For even n, the dominant contribution comes from the graph made of injl — 1) 
four-point vertices, which reads 



where R is assumed to be much smaller than M . For odd n, the graph with one three- 
point vertex and (n — 3)/2 four-point vertices makes dominant contribution to W\J\ and 
its magnitude can be estimated similarly. As a result we find, for the case a in J takes the 
value am(-R) defined in (p^), the graphs with the lowest order in the coupling constants 
make dominant contribution to provided that the coarse graining scale satisfies 

iV < 10. If the dominant contribution of the a— integral comes from a narrower region 
than (p2D, the constraint on N can be even weaker. 

Thus the PDF of the coarse-grained field up to the first order of the couphng 
constants is given by 




g-(fci+-+A;„_i)2fi2 



(fci + fc2 + k^Y + M2 ■ ■ ■ (fci + ■ ■ ■ + fc„_i)2 + M2 




(39) 
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exp [-'I. {(^,)% + ^ ((^h)')) co= + |! ((,.„)'))(40) 

The PDF with R = P/{2tt) is depicted in Fig. |^. The probabihty ipji > ip^^f = 
O.llTTc is 43.8% and the root-mean-square (RMS) is 0.76Tc, which is much larger than 
the inflection point ip-mf and the maximum point (f- = 0.277Tc of the potential. This 
result suggests that there is non-negligible amount of soaking into the asymmetric phase. 
Therefore we can conclude that the electroweak phase transition does not proceed by the 
standard bubble nucleation picture. 

For comparison, we show the probability distribution function for the field coarse- 
grained over the inverse mass scale R = 1/M in Fig. |. In this case, since R is large 
enough, from the physical argument, we can infer that the fluctuation becomes small and 
approaches the Gaussian configuration. Indeed, the probabihty of > is extremely 
small and the RMS is 0.028Tc, much smaller than the inflection point, which is consistent 
with the above observation and 0. Table 1 shows how the probability (fn > ip-mf and the 

root-mean-square ^ {^'r) change as a function of the coarse-graining scale R. Roughly 
speaking, if i? > 40/5/(27r), we can neglect the role of these thermal fluctuations. In the 
previous work |^, the critical scale is about Rl — 80/5/ (27r). The relation between the 
coarse-graining scale in this paper and that in the previous paper is not so clear because 
the former is defined in the spherical Gaussian window function and the latter corresponds 
to that defined in the cubic top-hat window function. But, if the correspondence is decided 
by the equality of the volume, that is, 

rf3^__/(0; R) , (41) 



we find Rl ~ 3.0R. Therefore the classical lattice simulation of is in reasonable 
agreement with the present analysis. 



4 Concluding remark 

In this paper, we formulate the probability distribution function for the coarse-grained 
fluctuations, which can be evaluated by calculating perturbatively the connected gener- 
ating functional W^[J]. In application, we apply our formalism to the electroweak phase 
transition. We evaluate the probability distribution function of the field coarse-grained 
over 2ti / f3c so that we find that the thermal fluctuations called subcritical bubbles play a 
significant role and that the universe is not homogeneously in the symmetric phase at the 
onset of the electroweak phase transition. We also calculate the probability P^Pr > v^inf] 
and the root-mean-square at several coarse-graining scale. The result is depicted in Table, 
m. We can observe that as the coarse-graining scale is larger, the probability P^Pr > v^inf] 
becomes smaller and also the RMS smaller. Thus, the coarse-graining scale is very im- 
portant in investigating the property of thermal fluctuations. Unless the coarse-graining 
scale is specified properly, the discussion becomes meaningless and brings out confusion. 
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Table 1: Property of thermal fluctuations at several scales 



R 




P 






/5/(27r) ~ 0.159/3 




44% 


0.76Te 


10/3/(27r) 


~ 1.59/5 




22% 


0.15Tc 


20/?/(27r) 


~ 3.18/3 




12% 


9.8 X lO-^Tc 


30/5/(27r) 


~ 4.77/3 




5.9% 


7.3 X lO-^Tc 


40/3/(27r) 


~ 6.37/3 




2.6% 


5.8 X lO-^Tc 


50/5/(27r) 


~ 7.96/3 




0.99% 


4.8 X lO-^Tc 


60/3/(27r) 


~ 9.55/3 




0.32% 


4.1 X lO-^Tc 


70/?/(27r) 


~ 11.1/3 


8.4 X 10-2% 


3.5 X lO-^Tc 


80/3/(27r) 


~ 12.7/3 


1. 


8 X 10-2% 


3.0 X lO-^Te 


1/M 


~ 13.7/3 


1.6 X 10-=^% 


2.8 X lO-^Te 
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Figure 1: One-loop improved finite temperature effective potential of the Higgs field, 
VewW\- 
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Figure 2: Contribution to W[J] up to 0{iJ,,X). 
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Figure 3: Probability distribution function for the field coarse-grained over /3/(27r). The 
inflection point, (fir^i in the potential corresponds to 0.117/3 and the local maximum, 
0.277/3. 
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ure 4: That over the inverse mass scale. 
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